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Inspired by the model of Jarvinen and Kiritsis, we present a simple holographic model for the on set 
of chiral symmetry breaking at the edge of the conformal window in QCD in the Veneziano limit. 
Our most naive model enforces the QCD two loop running coupling on a D3/D7 holographic brane 
system. The mass of the holographic field, describing the chiral condensate in the model, is driven 
below the BF bound when the running is sufficiently strong, triggering chiral symmetry breaking 
for Nf /N c < 2.9. This model though contains too great a remnant of supersymmetry and does not 
correctly encode the perturbative anomalous dimensions of QCD. In a second model we impose the 
QCD anomalous dimension result and find chiral symmetry breaking sets in for Nf/N c = 4 at a 
BKT-type phase transition. In this case the transition is triggered when the anomalous dimension 
of the mass operator 7 m = 1. 



I. INTRODUCTION 

There has been much interest in how the phase of QCD 
depends on the number of quark flavours for many years 
now. In the Veneziano limit, where the number of colours 
N c — > oo with fixed x = Nf /N c , we may treat a; as a con- 
tinuous variable. At x = 11/2 the one loop beta function 
vanishes. Just below that value of x the theory is known 
to be asymptotically free and to have a Banks-Zak fixed 
point [1, 2] at which the one and two loop beta functions 
balance to give a non-trivial, perturbative, conformal, IR 
fixed point. As x falls the fixed point value rises until the 
perturbative regime is lost. Based on the observation of 
chiral symmetry breaking in N c = 3, Nf = 3 QCD it is 
presumed that at some critical x c the IR conformal the- 
ory is replaced by one with a chiral condensate and a 
mass gap. 

A number of methods have been used to estimate 
x c . Truncated Schwinger-Dyson equations suggest 3.5 < 
x c < 4 [3, 4]. In these models chiral symmetry breaking 
is triggered when the anomalous dimension of the quark 
anti-quark operator hits of order one 1 [5]. The precise 
value for x c then depends on the truncation scheme, the 
choice made for the running coupling profile with energy 
scale, /x, and the anomalous dimension relation for the 
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1 In the Schwingcr Dyson analysis the criticality condition for 
chiral symmetry breaking is when 7 m (2 — "f m ) = 1 and hence 
j m = 1. To predict the critical value x c one needs to make an 
assumption for the dependence of 7 m on x at strong coupling. 
In [4] the one loop form is used and the criticality condition 
expanded at small 7^ giving 7^' = 1/2 as the condition if 
extended to strong coupling. Using the two loop form for the 
running coupling in 7^ gives the prediction x c = 4. 



quark mass term, j m . Other attempts to estimate the 
critical value have been made in [6]- [11] and typically give 
a similar estimate. Recently there has been much inter- 
est in simulating such theories on the lattice too [12]-[21] 
the simulations are still in an early phase but already 
support the general picture from all these analyses. 

Jarvinen and Kiritsis [22] have recently proposed a 
holographic model of the strongly coupled near con- 
formal regime around x c (the work in [23-28] is also 
very relevant). Their model consists of 5D supcrgrav- 
ity with a dilaton field dual to the running coupling and 
a "tachyon" field dual to the chiral condensate. They 
impose potentials for all these fields that generate the 
known two loop running for the coupling and the pertur- 
bative relation for the anomalous dimension in the UV of 
their description. They predict the range 3.7 < x c < 4.2 
and that the transition at Xq IS cL BKT type transition 
in which the condensate grows exponentially from the 
transition (as expected - see [29, 30]). In a sense all 
this physics is imposed by the choice of potentials in the 
model but those choices are reasonable and it is encour- 
aging that the results match other estimates. 

In this paper we wish to attempt a similar construction 
using an alternative holographic model of chiral symme- 
try breaking. The D3-D7 system [31] provides the sim- 
plest and best understood holographic description of a 
strongly coupled gauge theory with quark fields. At heart 
it consists of SU(N C ) AT = 4 super Yang-Mills theory 
with Nf M = 2 quark hypcrmultiplcts. In the quenched 
approximation the theory is conformal and on the grav- 
ity side is described by probe D7 branes in A0IS5 x S 5 . 
The theory is 3+1 dimensional at all energy scales and 
has a conformal UV in which the identification of the op- 
erator matching between the field theory and the gravity 
description is clean. The simplest description of chiral 
symmetry breaking is found by imposing a background 
magnetic field on this theory [32] - the description is reg- 
ular throughout and the interpretation again clear cut. 
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Interestingly the DBI action for the probe with the mag- 
netic field present is equivalent to the same theory with 
a particular choice of running gauge coupling. This effec- 
tive dilaton is not backreacted on the geometry. It there- 
fore seems natural to move to phenomenological models 
where one simply imposes some running coupling on the 
theory by hand - the underlying reaction of the holo- 
graphic description seems likely to correctly capture the 
resulting physics of chiral symmetry breaking. Indeed 
in a recent paper we have looked at the phase struc- 
ture of just such a model with a running coupling with 
a step change between two conformal regimes [33]. The 
imposition of the two loop QCD running is a very simi- 
lar analysis which we explore here, concentrating though 
on the transition to chiral symmetry breaking. Placing 
the probe brane in the presence of the dilaton matching 
the two loop gauge running essentially looks at the dy- 
namics of one of the Nf quark flavours in a background 
backreacted to the full dynamics of the Nf quarks. Our 
model is more direct than that of Jarvinen and Kiritsis 
[22] in that we simply impose the running of the gauge 
coupling, and also in a later analysis the QCD anomalous 
dimension relation, rather than imposing a potential and 
then solving for these quantities. If one had the correct 
gravity dual of the gauge theory then the more involved 
process would capture more of the dynamics but if we 
are simply modeling the gauge theory then our approach 
may be sufficient. 

First we will very naively impose the two loop QCD 
beta- function on the D3/D7 dynamics. We find that chi- 
ral symmetry breaking is induced for x < 2.9 and that 
the transition to the chiral symmetry breaking phase is 
second order in nature. This value of x c is low relative to 
other estimates and the transition type is at odds with 
that argued for in QCD in [29, 30]. To understand this 
we recast the DBI action for small fluctuations about the 
chirally symmetric phase as a slipping mode in AdS$ (we 
study this analysis in more detail in the Appendix) . One 
can then plot its mass squared as a function of the ra- 
dial coordinate and seek a violation of the Breitenlohner 
Freedman (BF) bound [34] which would lead to an insta- 
bility. The model only shows an anomalous dimension 
for the chiral condensate in the regimes where the dila- 
ton is running and the size of the anomalous dimension 
is proportional to the strength of that running. Our crit- 
ical value of x c = 2.9 corresponds to the theory which 
first has sufficiently strong running present. In terms of 
the slipping mode mass squared the BF bound must be 
violated over a sufficient interval in the radial direction 
of the gravity description - within that interval the BF 
bound is substantially violated at the transition leading 
to the second order behaviour. 

This analysis highlights a failure of the D3/D7 sys- 
tem as analyzed so far - it has too much supersymme- 
try present. In the IR conformal regime the background 
gauge dynamics returns to that of N = 4 Super Yang- 
Mills. It has too much symmetry and does not induce an 
anomalous dimension for the quark mass / condensate no 



matter how large the gauge coupling. This is in direct 
contradiction to QCD where the anomalous dimension 
7 m is directly proportional to the magnitude of the cou- 
pling, at least in the perturbative regime [35]. Simply 
put we need to introduce more supersymmetry break- 
ing into the description. We show how by a choice of 
background dilaton the QCD one-loop anomalous dimen- 
sion relation can be imposed on the model by hand. We 
next impose on top the two loop QCD running profile 
within the anomalous dimension relation. In this model 
the slipping mode's mass squared, m 2 , asymptotes to —3 
in the UV and to some lower IR fixed point value. As 
it passes through the BF bound of —4 chiral symmetry 
breaking is triggered. The two loop running's IR fixed 
point implies that at the transition the IR mass squared 
lies at exactly —4 and this is the condition needed for 
a BKT transition (see [36, 37] for the first examples of 
holographic BKT transitions), which we indeed observe. 
In this model x c = 4. 

It's worth stressing that the BF bound is violated in 
these holographic models precisely when m 2 = —4 and, 
using the usual conformal AdS mass-operator dimension 
relation, 7, n = I. This seems a robust holographic predic- 
tion, particularly since we envisage a conformal IR regime 
where the AdS mass-dimension relation is expected to 
hold. Note that x c and the BKT transition behaviour are 
completely determined by the IR fixed point behaviour of 
the coupling and the precise non-perturbative running is 
not crucial. There are more vagueries in the precise pre- 
diction of x c since we must assume a non-perturbative 
relation between the anomalous dimension and the value 
of the IR coupling. We have used the leading perturba- 
tive relation between m 2 and the one loop anomalous di- 
mension and extended it to the non-perturbative regime, 
giving criticality when 7m =1/2 and x c = 4. Given the 
full QCD dynamics this value may be different though. 



II. THE D3/D7 SYSTEM 

Our starting point is the holographic D3/D7 system 
[31]. Strings tied to the surface of the iV c D3 branes 
generate the adjoint representation fields of the N~ = 4 
gauge theory. Strings stretched from the D3 to the D7 
are the quark fields in the fundamental representation of 
the SU(N) group. 

In the strong coupling limit the D3 branes are replaced 
by the geometry that they induce. We will consider a 
gauge theory with a holographic dual described by the 
Einstein frame geometry AdSs x S 5 




where we have split the coordinates into the X3+1 of the 
gauge theory, the p and which will be on the D7 brane 
world-volume and two directions transverse to the D7, 
w^jWq. The radial coordinate, r 2 = p 2 + w 2 + Wg, cor- 
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responds to the energy scale of the gauge theory. The 
radius of curvature is given by i? 4 = 4ng s Na 2 with N 
the number of colours. The r — > oo limit of this the- 
ory is dual to the Af = 4 super Yang-Mills theory where 
9s = .9uv ^ s t ne constant large r asymptotic value of the 
gauge coupling. 

In addition we will allow ourselves to choose the pro- 
file of the dilaton as 1 — > 0. Simplistically this represents 
the running of the gauge theory coupling, = /3, where 
the function /3 — > 1 as r — ¥ oo. For the coupling pro- 
files we will consider later the UV form of j3 will have 
weak logarithmic running present - we will impose a UV 
cut off when f3 = 1 corresponding roughly to the scale 
where the holographic dual should be matched to per- 
turbativc QCD. Above that cut off we simply set /3 = 1. 
The physics we study is all in the IR and not affected 
by the precise form of this cut off though. In the final 
section we will use the dilaton function as an input to 
the DBI action for the quark physics to enforce the QCD 
anomalous dimension relation. At this point the rela- 
tion between the dilaton, the gauge coupling and phe- 
nomenological corrections to the DBI action become less 
clear but our philosophy is simply to phcnomcnologically 
enforce the correct quark physics in the DBI action. 

We will introduce a single D7 brane probe into the ge- 
ometry to represent the dynamics of one quark in the 
theory - by treating the D7 as a probe we are working in 
a quenched approximation although we can reintroduce 
some aspects of the N f quark loops through the running 
coupling's form. This system has a U(l) axial symmetry 
on the quarks, corresponding to rotations in the w 5 -wq 
plane, which will be broken by the formation of a quark 
condensate. The hope is that the dynamics of chiral sym- 
metry breaking for the quark described by the probe is 
generic across many gauge theories and the results will 
be applicable to QCD. 

We find the D7 embedding function e.g. 11)5(0), wq = 0. 
The Dirac Born Infold (DBI) action in Einstein frame is 
given by 

S D7 = -T 7 J d 8 £eV-detP[G] o6 

= -T 7 J d 4 x dp p 3 ^! + {d p LY , 



(2) 



where w 5 = L, T 7 = (2tt)- 7 a'- i g^ and T 7 = 2tc 2 T 7 
when we have integrated over the 3-sphere on the D7. 
The equation of motion for the embedding function is 
therefore 



Pp 3 d p L 

V 1 + 



2V V / 1 + (W|§=0. (3) 



The UV asymptotics of this equation, provided the dila- 
ton returns to a constant so the UV dual is the N = 4 
super Yang-Mills theory, has solutions of the form u> 5 = 
d + c/p + • ■ ■ , where we can interpret d as the quark 
mass [m q = c?/27ra') and c is proportional to the quark 



condensate. 

The embedding equation (3) clearly has regular solu- 
tions u>5 = m when f3 is independent of r - the flat embed- 
dings of the Af = 2 Karch-Katz theory. Equally clearly 
if d/3/dr 2 is none trivial in then the second term in 
(3) will not vanish for a flat embedding. 

There is always a solution w§ = which corresponds to 
a masslcss quark with zero quark condensate (c = 0). In 
the pure Af = 2 gauge theory with j3 = 1 this is the true 
vacuum. In the symmetry breaking geometries [33, 38] 
this configuration is a local maximum of the potential. 

If the coupling is larger near the origin then the D7 
brane will be repelled from the origin ending at p = 
with L'(Q) = 0. The symmetry breaking of these solu- 
tions is visible directly [38]. The U(l) symmetry corre- 
sponds to rotations of the solution in the w^-wq plane. 
An embedding along the p axis corresponds to a mass- 
less quark with the symmetry unbroken. The symmetry 
breaking configurations though map onto the flat case at 
large p (the UV of the theory) but bend off axis breaking 
the symmetry in the IR. L(0), the IR quark mass, is a 
good order parameter for studying the chiral symmetry 
breaking that also reflects the bound state masses of the 
theory. 



III. IMPOSING THE 2-LOOP QCD RUNNING 

Our first analysis is straightforward. We impose the 
two loop running of the QCD gauge coupling on the dila- 
ton profile of the D3/D7 system. That running is deter- 
mined by 



where 



dX 
dp 



2 (II - 2x) 

3 (4^) 2 ' 



-b Q X 2 



hx 3 



2 (34 
3 



13z) 



(4tt) 4 



(4) 



(5) 



In the bi, we omitted a subleasing term 0{N~ 2 ) at large 
N c . We simply identify the radial direction r with the 
RG scale of the field theory and set A = j3. As is well 
known these equations have logarithmic running in the 
UV and an IR fixed point that grows from zero as x is 
reduced from x = 5.5 2 . 

Here the UV is not strictly conformal although it ap- 
proaches it at weak coupling asymptotically. Neverthe- 
less it is easy to look for chiral symmetry breaking. We 
continue to associate massless quarks with D7 embed- 
dings that approach the p axis at large p and seek so- 
lutions that bend off axis with that UV boundary con- 



2 Note for reference that in the usual gap equation analysis [3, 4] 
the critical coupling is given by \ c = 8n 2 /3 which is first achieved 
in the IR for x c = 4 
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(a)Thc model of section III where we naively 
impose the QCD running coupling and shows 
a second order transition at x = 2.949. 




(b)The model of section IVa where we impose 
the IR QCD anomalous dimension relation 
and shows a BKT transition at x c = 4. 

FIG. 1: Plots of the IR mass L(0) against x = N f /N c for 
our two models. 



(a)The model with the QCD running 
imposed in section III (x = 3.5, 3.29, 3.0). 
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(b)The model of section IVa where the QCD 
anomalous dimension is imposed in the IR 
(x = 4.5,4,3.5,3.3). 

FIG. 2: Plots of the AdSs slipping mode m against r in our 
two models. 
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dition. In fact the simplest identifier of chiral symmetry 
breaking is to look for solutions that begin with £'(0) = 
and shoot out to lie below L = in the UV. We use the 
value of L(0), the IR quark mass, as the order parameter 
for chiral symmetry breaking. 

In Fig la we show a plot of L(0) vs x. The transition 
is clearly second order and by fitting we determine it to 
be mean field with critical exponent 1/2. This second 
order nature is of course at odds with expectations that 
the transition with x should be of the BKT type [29, 30]. 

To understand this behaviour let us perform a lin- 
earized analysis on our DBI action to see why the flat 
embedding L = becomes unstable. 

We have an action, which is proportional to (2), 



S= / dpA(r)pVl + L" 2 



(0) 



where r 2 = L 2 + p 2 . We expand for small L 



Hp)p 3 - 



'1 1 



dp ^ dp ' 



Jp Ap3 



Now the first term in our action can be recast by set- 
ting L = p4> 

S = J &p l - P 3 L 2 = J dp 1 - (W 2 - 3p 3 2 ) , (9) 

where L' = dL(p)/dp. This is the action of a canonical 
in 2 = —3 scalar in AdSs. The remaining term in the 
action becomes 

s = [ dphift 2 . (io) 

J 2 p dp 

So we have a AdSs scalar with p dependent mass squared 



S 



dp ( 2 A ( r ) 



P 3 L' 2 +p 3 



L=0 



dX 
dL 2 



L- 



L=0 



(7) 



where L' = dL(p)/dp. To make the kinetic term canon- 
ical, we can now make a coordinate transformation 3 on 



= —3 — 5m 2 , Sm 2 — —A 



ffd\ 
J* dp 



(11) 



We plot this mass against p in Fig 2a for a variety of 
choices of x. We first note that the mass squared of the 



3 See Appendix for more detailed and general discussion on the 



coordinate transformation. 
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slipping mode approaches —3 in both the UV and IR, 
which we will return to shortly. The instability that is 
causing our phase transition with x is in the intermediate 
period where the mass squared is falling below the BF 
bound of —4. Note that x = 3.29 is the first case where 
the BF bound is met at one point in p but that this 
value of x is significantly above the critical value x c = 
2.95 found above. Apparently the BF bound must be 
violated in a region of p for the instability to trigger a 
transition. At the point of transition the BF bound is 
violated in a range of p and the effective mass squared 
lies considerably below the BF bound in the mid-region. 
Such a scenario has caused a second order transition. See 
(A15) and footnote 4 for more examples. 



IV. IMPOSING THE QCD ANOMALOUS 
DIMENSION 

Our naive model above of the x dependence of QCD 
suffered from an excess of supersymmetry in the IR 
regime, left over from our underlying construction. The 
model included the running gauge coupling but not the 
QCD anomalous dimension relation. We will now enforce 
the perturbative QCD form of that relation (with the two 
loop QCD form for A) on the model as an alternative way 
to include the QCD physics. That is we use the two loop 
relation to fix the coupling at the IR fixed point and then 
use the one loop anomalous dimension relation to predict 
x c from the point where the slipping mode mass becomes 
-4. 



This plot of the mass squared of the slipping mode 
against p reveals a number of failings of our most naive 
model. In particular the mass squared returns to —3 in 
the IR. The reason is that the gauge coupling becomes 
constant in the IR conformal regime and the brane con- 
struction returns to that of the N — 4 gauge theory with 
quarks. The mass squared is —3 because the model re- 
turns to a highly supcrsymmetric configuration in the IR 
where the anomalous dimension of qq is protected to be 3. 
This is quite unlike in non-supersymmetric QCD where 
the anomalous dimension, j m , of the quark mass m q (the 
dimension of qq is 3 — 7m ), at one loop, is given by 



(i) 



d In m„ 



/'- 



df i 



3A 
(4^ 



(12) 



At the IR fixed point one expects a non-zero j m . How- 
ever, in our holographic model, using the naive scalar 



mass operator dimension relation, 
find 



nr 



A(A - 4), we 



A. IR Physics 

Let us first consider the IR conformal regime where 
we want a constant non-zero value of 7 /_r. Our model 
predicts the slipping mode mass (11) 



L IR 



p 5 dp 
p 4 dp 



(15) 



where we substitute /3 for A since we have given up the 
identification of /? in the DBI action with the gauge cou- 
pling. Here we are concentrating on making j m match 
QCD instead so the quark physics is correct. The choice 
of /3 that gives such a constant m 2 IR is (3 
for which we find 



js (0 < q < 2) 



m IR = — 3 — Sm 2 



Sm 2 = 



(2-g) 2 



(16) 



By (13), Sm 2 is related to j„ 



7m = 3 - A = 1 - v 7 ! - Sm 2 (13) 
Sm 2 = l-(l- 7m ) 2 , (14) 

where Sm 2 is defined in (11). Strictly the relation m 2 = 
A(A — 4) is valid only in conformal regimes where the 
scalar mass is constant but we allow ourselves to use it 
slightly more liberally here. Therefore, the model we 
present only conjures an anomalous dimension in the 
regime in which the coupling is running (Sm 2 ^ 0), 
breaking both conformal invariance and supersymmetry. 
For a model that is precociously asymptotically free such 
as x = 1 QCD this deviation of our model from QCD 
is probably not so important for the phenomenology - in 
both cases the coupling grows rapidly and a quark con- 
densate is triggered. If we wish to model the transition to 
chiral symmetry breaking though with changing x, where 
we leave an IR fixed point theory, it is more crucial. Our 
estimate of x c = 2.95 is most likely an under estimate be- 
cause we have not included the contribution to "f rn from 
the absolute value of the coupling A. 



1-4 1- 



[2-qf 



(17) 



Here the use of (13) is more valid than in the previous 
section; the scalar mass in the action given by (9) and 
(10) is constant and by ansatz p matches the RG scale 
of QCD. It is possible that back-reaction between the 
geometry and the scalar might disturb this relation but 
it seems fairly sound. Note the conditions that m 2 R = —4 
and 7 m = 1 are the same,where q = 0.536. 

This model displays a BKT transition as q is changed 
continuously through 7m = 1. As usual the BKT transi- 
tion occurs due to the presence of an infinite number of 
unstable, Efimov modes at the transition point [30]. We 
can see them here explicitly by considering the action for 
static, linearized, mesonic solutions around the L = 
embedding. The action is 



S = 



F 2 



PL 



'2 



4^ 2 
p 4 
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where a prime is a p derivative and dot a time derivative. 
If we move to the inverse z-coordinate (z = 1/p) we have 



result 



S 



y -5 



dz— [ ?z 4 L 2 



Pz *L 2 - z 4 ^L 2 
oz z 



(19) 



where a prime is now a z derivative. We can now write 

the equation of motion for a solution of the form L = 

e -iut z (i- q )/2^ and p = z q 

(3-8g + g 2 ) 1 2 

- ip + = u ip , 20) 

4 z z 

which is a ID Schrodinger equation form with a 1/z 2 
potential. This problem is known and becomes unsta- 
ble when the coefficient of the 1/z 2 term is equal to 
— 1/4. This condition is equivalent to j m = 1 in (17). 
At that point an infinite number of unstable negative en- 
ergy modes emerge from E = 0. At the critical value 
of — 1/4 all of those modes play a role in the transition 
generating the BKT transition. 

This discussion so far has been restricted to the IR and 
a more complete model would require that /? — > 1 in the 
UV. A simple fix is to set j3 = 1 + c/p q . In this case, 



6m 2 (p; qi c) = cqp- 2 i(c + pi) 2 F? 

Its IR asymptotic behavior is 

4g ( _ 1 



2 2 
1,- - 

q 



-cp 



5m' 



(2-9) s 



1 - 



o(i -qV 



(21) 



(22) 



which is the same as (15) with a p-dependent correction. 
The IR behaviour matches our discussion above. 

At this point we can make a simple model to extract 
the critical value of Nf in QCD. The two loop QCD beta 
function has a fixed point at 



A,„ 



11 -2x 
13x - 34 



(4tt) 2 



In the Banks-Zak regime where perturbation theory ap- 
plies, 7 m * = . the order A relation between the Sm 2 
and 7 m * is given by 



2-v (1) - 



6A» 
(4tt) 2 



(24) 



where we used (14) and (12) and 7m* denote the order A 
relation. 

Of course we have no true idea how to continue this 
relation into the non-perturbative regime but following 
the spirit of [4] we will simply assume (24) applies at all 
values of the coupling. The holographic model tells us 
that the transition will occur when m 2 = —4 (8m 2 = 1) 
so we find, using the one loop QCD anomalous dimension 



6A* 

W) 



, 11 - 2x 
' 13a;- 34 



(25) 



This gives x c = 4. Note that this amounts to 7m* = 1/2, 
which coincides to the one-loop pcrturbative field theory 
computation [4]. 

Finally we can numerically check the BKT nature of 
the transition as well. We can simply set A = 1/p 9 with 
q and x related, through the IR relations (16) and (24), 



4q „ 11 - 2x 
6- 



(2-<z) 



13x - 34 



(26) 



We then numerically solve for the D7 embedding, L as 
a function of p. L(0), the IR quark mass, is a useful 
order parameter - we show the result for L(0) vs x in Fig 
lb - the BKT type transition is apparent with x c = 4. 
Close to x c this simple model and the case j3 = 1 + 1/p 9 
coincide since the dynamics is dominated in the far IR. 



B. All RG scales 

To construct a full model of the RG flow in the con- 
formal window, one should enforce the QCD anomalous 
dimension formula (24) at all energy scales or p. In par- 
ticular we want 



p 4 dp 



6A(p) 
(4tt) 2 



(27) 



To find the associated /3 one can re-arrange for p, 




(28) 



(23) differentiate, and find the differential equation 



2 



-6A 



(4tt) 2 p 5 /3^' 



1/2 



(29) 



We can solve for (3 numerically by shooting from some 
initial value of p and trialling various values of the ini- 
tial condition j3'. Typically the true solution lies on the 
crossover between solutions that are real at all p and 
those that go complex so the correct initial condition can 
be tuned to. Once found the numerical solution can be 
tested that it is a good solution of (29) and that it has the 
IR fixed point behaviour 1/p 9 where q and x are related 
by (26). 

We can then use these solutions to solve for L as a 
function of p - close to x c the results are again those in 
Fig lb since the dynamics is entirely determined by the 
IR fixed point. 
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V. SUMMARY 

We have presented two simple holographic models of 
x = Nf/N c behaviour of QCD at large N c . In our first 
model we imposed the QCD two loop running directly 
on the D3/D7 system through a non-backreacted dilaton 
profile. We found chiral symmetry breaking sets in at 
x c = 2.95 at a second order transition. The transition is 
expected to be at a larger value of x and to be of BKT 
type [29, 30] and we highlighted that this discrepancy 
is due to the IR supersymmctry of the model forcing 
7, fl = 0. In a second model we imposed the perturbative 
QCD 7,„ relation and found a BKT transition at x c = 4. 

Whilst these models are much less sophisticated than 
the very nice model of Jarvinen and Kiritsis [22] , in which 
the AdS-space backreacts to the running coupling and 
the quark condensate, we believe they highlight the key 
ingredients. One must input into the model, either di- 
rectly as we do, or indirectly through supergravity po- 
tentials as in [22] , the form of the running coupling and 
the impact that has on the quark anomalous dimension. 
Since we do not have the true QCD dual all of this is the 
model builder's choice. The clear prediction from AdS is 
that the chiral transition will occur when the AdS slip- 
ping mode associated to the quark condensate hits a mass 
squared at the IR fixed point of —4, the BF bound. This 
corresponds to j m — 1. The Miransky scaling or BKT 
nature of the transition is then also very clear in the 
holographic description through the presence of Effimov 
modes. 



where r(p) = y/ L 2 + p 2 and d is an integer related to 
the dimension of the background and worldvolume space- 
time. For example, for the D7(D5) probe brane in AdSs 
x 5 5 , d = 4(3). We assume that 



(q<d-l) p 



OQ 





CUV) 
(IR), 



(A2) 



where c is constant. The first condition comes from the 
fact that the slipping mode (0 = L/p) is a scalar in 
AdSa+i in UV. The second condition restricts us to an 
effective IR AdS space. When q = d — 1, IR space is 
effectively AdS?- 

The linearized action in terms of the slipping mode 
reads 



d - 1 (pV 2 

m = (1 - d) - — 

dlogp 

where O = P(r(p), p)\ L =o- 



m 2 (j> 2 ) 



d(3 
A) dL 2 



(A3) 



L=0 



1. Effective geometry changed 

In the UV, (3q = 1 and the action (A3) corresponds to 
the scalar action in AdS^+i space with the UV mass 



l uv 



l-d> -tf 2 /4, p -> oo . 



(A4) 
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For all d, the BF bound is satisfied (in AdSd+i) 
In the IR, the action is written as 



S 



,2j2 



(l-d + q) + 2 



p 2 dp_ 
Jv dL 2 " 



L=0 



(A5) 
(A6) 



The scalar effectively lives in AdSd- q , where q < d — 1 . 
To go further we consider two cases: j3 = /3(r(p)) and 
p = P(r{p),p). 
case 1: For (3 = (3(r(p)), 



Appendix A: Effective scalar mass and the BF mass 
violation 



df3_ 
dL 2 



L =o 2 P d P ' 



(A7) 



In this appendix, we show how to identify the effective 
mass of the slipping mode of the probe brane in an ef- 
fective AdS space in a more general context and in more 
detail. 

In general the action of the embedding L(p), which 
is a function of only p, a holographic direction, can be 
written as 



S= / dpP{r{p\p)p d - 1 ^l + U{p) 



(Al) 



L IR 



(l-d + q) + 



d\og/3 
d\ogp 



1 



(A8) 



Note that the m 2 is the same in UV and IR. However, the 

stability criteria, the BF bound — ^ d ~^ , is now changing 
and a function of q. Therefore, if q is a continuous param- 
eter (for this purpose, let us continue q to real values), 
then the BKT transition occurs at q = d — y Md — 1). 
For d = 4, q ~ 0.536, which is the same value we ob- 
tained in section IV. 
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case 2: For (3 = /3(r(p), p), we have to study case by case, 
since (A7) is not valid. As an example, let us consider 
D3/D7(D5) at finite B and density, d [36, 37, 39, 40]. 



be of mean- field type 4 . (The infinitesimal violation of 
the BF bound as in (All) is a characteristic of the BKT 
transition.) 



a. D3/D5 



S = / dpP(r(p),P)p 2 Vl + L>( P y 



where 



P(r(j>),p) 




(A9) 



(A10) 



By (A4), the UV mass is -2 in AdS 4 , while, by (A6), 
the IR mass is (d — 3, q = 2) 



2. Effective geometry fixed 

There is alternative way, in which we keep the UV 
AdS space for all p. For this, wc need to redefine the 
coordinate system by 



so 



dp 
dp 



/V" 



pd-i 



(A16) 



2B 2 



l IR 



d 2 + B 2 



(All) 



in AdS2- The BF bound is violated at d = \flB and this 
violation by the continuous parameter, d or B, implies 
the BKT transition. If d = then mf R = —2, which is 
consistent with (A8). 



b. D3/D7 



S = / dp(3(r{p),p)p^l + U(p) 2 . 



where 



P(r(p),P) = 




(A12) 



(A13) 



By (A4), the UV mass is -3 in AdS 5 , while, by (A6), 
the IR mass is (d = 4, q = 3) , 



B 2 p 2 
d 2 



l IR 



0. 



(A14) 



in AdS2 ■ It satisfies the BF bound for all B and d. 
However, the instability is in the intermediate regime. 
We can see this by expanding the action in the regime 
d/5<p<C \/~B, of which linearized equation of motion 



L" + -L' + 2^L = 0. 
P P 



(A15) 



The slipping mode is effectively the scalar of m 2 = — 3 in 
AdS3, which violates the BF bound. It also can be seen 
more directly from (A3), where the second term and third 
term cancel out, leaving the first term, (1— d) = —3. Note 
that this instability happens only for a large enough B (or 
small enough d) to satisfy the condition d/B <C p vB, 
Note also that the BF mass violation is finite as we dial 
B for a fixed d, and the phase transition turns out to 



(A17) 



which is defined only for < q < d — 2. Otherwise, the 
integral diverges and p is not defined. (So, the previous 
examples of the D7(D3) probe in AdSs x S 5 cannot be 
analyzed in this way; q is too big.) In terms of a new 
coordinate p, the action reads 



S 



dp 1 -^ 1 (pV 2 



mV) 



m 2 = (1 - d) + 



J2d-3 



(A18) 



P 



2d-4 



where we consider only the case /3 = /3(r(p)), so that we 
can use (A7). Note that for a function /3q = p~ q (the 
normalization of j3o does not matter, since any normal- 
ization factor is canceled in (A18)), rri 1 is constant: 



m 2 = (1 - d) - q 



d-2 



d-2 



(I 



(A19) 



One might wonder if this analysis is consistent with 
the previous one (Appendix A.l). For example, for d = 
i,q = 1, both analyses are applicable. They must be 
consistent since the BF bound analysis is an effective 
tool and how to interpret the action should not change 
the physics, i.e. for d = 4, q = 1, we can interpret the 
action of cither (1) a scalar in AdSs with m 2 = —7 or (2) 
a scalar in AdS4 with m 2 = —3. However, both cases tell 
us the scalar mass violate their own BF bound, so they 
are consistent. To see this more clearly, let us check the 



4 A non-mean field (but non-BKT) type transition also can be 
understood in the same way. In the model studied in [41], 
it can be shown that an instability can arise in the range 
({(S/0) 1 /< 2 - A ),(d/0) 1 /( 2 - A )} « p < O 1 ^), where O is a 
phenomenological operator with dimension A. This range is es- 
sentially where the operator O dominates over B, d. 
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BF bound conditions, which arc 

l-d=_£^£ (A20) 
in Appendix A.l, and 

in this subsection. These seemingly different conditions 
indeed give the same results: the BF bound is violated 
at the value of q = q c , 

q c = d- 2Vd - 1 . (A22) 

Therefore we may interpret our analysis as either (1) 
m 2 does not change but the effective background is 
changing (2) m 2 is changing but the geometry does not 
change. 

Of course, we can do a mixture: partial change of ge- 
ometry and partial change of m 2 . How this works in gen- 



eral can be seen by the following simple example. The 
equation of the scalar field at the boundary of AdS^+i 
space (z 0) reads 

$" + - -$'-m 2 — =0. A23) 

z z z 

By the definition $ = z^~(j>, it can be transformed to 

,/ (l-<0 „ ( 9 d 2 -D 2 \ 4> . 
^ + L V J ^' - (m 2 - _j *=0. (A24) 

It is formally interpreted as the scalar in AdS<j+i space 
with the modified mass m 2 — (d 2 + D 2 )/4. In both cases 
the BF bound is the same, m 2 = — D 2 /A, so the physics 
does not change. Especially, the equation for D = 2 
(or d = 2) corresponding to the effective AdS2 is the 
Schrodinger equation with the 1/z 2 potential term and 
the — l/4z 2 potential plays a role for the BKT transition 
as discussed in [30] and section IV. 
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